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ABSTRACT 

We present new Spitzer photometry of the Eagle Nebula (M16, containing the optical cluster NGC 
6611) combined with near- infrared photometry from 2MASS. We use dust radiative transfer models, 
mid-infrared and near-infrared color-color analysis, and mid-infrared spectral indices to analyze point 
source spectral energy distributions, select candidate young stellar objects (YSOs), and constrain their 
mass and evolutionary state. Comparison of the different protostellar selection methods shows that 
mid-infrared methods are consistent, but as has been known for some time, near-infrared-only analysis 
misses some young objects. We reveal more than 400 protostellar candidates, including one massive 
young stellar object (YSO) that has not been previously highlighted. The YSO distribution supports 
a picture of distributed low-level star formation, with no strong evidence of triggered star formation 
in the "pillars". We confirm the youth of NGC 6611 by a large fraction of infrared-excess sources, 
and reveal a younger cluster of YSOs in the nearby molecular cloud. Analysis of the YSO clustering 
properties shows a possible imprint of the molecular cloud's Jeans length. Multiwavelength mid-IR 
imaging thus allows us to analyze the protostellar population, to measure the dust temperature and 
column density, and to relate these in a consistent picture of star formation in M16. 

Subject headings: infrared: stars - methods:data analysis - stars:formation 



1. INTRODUCTION 

The fusion energy generated in stars dominates the 
evolution of galaxies, heating and processing their inter- 
stellar media. However, our understanding of how stars 
form is incomplete. Some fairly reliable empirical scal- 
ing laws exist, and a believable scenario exists for how 
an isolated low-mass star can form via disk accretion, 
but we lack a detailed understanding of how an entire 
molecular cloud turns into a cluster of stars. Most stars 
form in cluster environments, and stellar densities imply 
some interaction between YSOs during their formation. 
We must therefore analyze entire star forming regions, at 
long enough wavelengths to not be limited by extinction, 
and determine the physical properties of whole popula- 
tions of young stellar objects (YSOs). Trends in evolu- 
tionary state, mass, and spatial distribution can help dis- 
tinguish between star formation scenarios such as com- 
petitive accretion (e.g. iBonnell et aLll2004D and predes- 
tination (in which cores determine stellar masses, e.g. 
iMcKee fc Ta^[200l . 

One particularly vexing aspect of the complex and 
important topic of star formation is that of triggering. 
The formation of massive stars feeds energy back into 
the nearby interstellar medium, compressing, irradiat- 
ing, and heating the natal molecular cloud. This feed- 
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back can have destructive or constructive effects, but it 
is not clear which dominates in a given cloud or overall 
in a galaxy. Although many examples exist of regions 
in our Galaxy and other galaxies with suggestive spatial 
distributions of young and somewhat younger stars, de- 
termining the relevant physical conditions and timescales 
to make a strong argument for triggering is more difficult. 
For some examples of these argumen ts, see the discus- 
sion of the Gal actic region W3/ 4 bvlOev et all ([2005), 
30 Doradus by IWalborn et al.l (120021)7 " various regions 
and theory by iDeharveng et al.l (I2005D. the p articular 
case of cloud-crushing described bv lLefloch et a l. ( 1997}, 
and further theoretical considerations in Elme green et aLl 
(j2.002). This issue must be clarified by making detailed 
assessment of whether triggered star formation is hap- 
pening in particular molecular clouds. 

The Eagle Nebula (M16, containing the optical clus- 
ter NGC 6611) contains a well-studied region of pos- 
sible triggered star formation, in the widely publicized 
"pillars of creation" o r "ele phant trunks" imaged with 
HST bv lHester etaH ([19961 hereafter H96). NGC 6611 
is a yo ung open cluster, l ikely still in its formation 
stages . iHillenbrand et al.l ()1993D and Ide Winter et al.l 
(i l997i ) found a spread in ages of the optically visible clus- 
ter members: The massive (M^, > IOMq) stellar popula- 
tion is ~2±1 Myr old, with the most massive members 
apparently having begun to evolve off of the ZAMS, and 
evidence o f at least one evolved (6 Myr old) 3OM0 star 
((Hillenbra nd et aT1ll993f ). There are also optically vis- 
ible (no longer embedded) intermediate mass stars (3- 
8Mq) that are still evolving onto the ZAMS. The fa- 
mous pillars are dense knots of molecular gas (jPoundl 
|1998[) and dust which shield the region behind them from 
the destructive radiation of the O stars in the central 
cluster. H96 identified "Evaporating Gaseous Globules" 
(EGGs) in and around the pillars, and a central ques- 
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tion has been whether the EGGs are being compressed 
and forming stars fast enough before they are photo- 
destroyed. Detailed near-infrared (NIR) searches find 
that about 15% of the EGGs have evidence of YSOs 
jMcCaughrcan fc Andersen 20021: [Thompson et al.ll2002t 
ISugitani et al.. 2002). Millimet er continuurn and rnolecu- 
lar lines (|Andersen et ani2004l : IWhite et al.lll99llPoundl 
119981) reveal only one or two possibly younger objects in 
the pillars and at their base. 

Mid IR imaging offers new insight into this region 
in at least two major ways. First, not all objects 
with protostellar dust show near-infrared (JHK) excess 
(e.g. only half of YSOs w ith 10 ^m excess in Taurus - 
Auriga have NIR excess; iKenvon fc HartmannI I1995D . 
and modeling protostellar spectral energy distributions 
(SEDs) is much more reliable with longer wave length 
data (e.g. lWhitnev et al.ll2004bl : lAllen et al.ll2004 ). Sec- 
ond, the unprecedented mapping speed of Spitzer al- 
lows imaging of a much larger are a than the previ- 
ous NICMOS ( Thompson et al."'2GG2 l) and ground-based 
(jSugitani et al.. 2002; McCaughrca n fc AndersenI l2002f) 
studies centered on the pillars, permitting an investi- 
gation of the protostellar population in the entire re- 
gion. This comprehensive assessment reveals star for- 
mation throughout the region, with the least evolved ob- 
jects most highly clustered together. We also use mid-IR 
data to map the dust temperature and average extinc- 
tion; the YSOs are preferentially located in regions of ele- 
vated dust temperature, and aside from NGC6611 which 
has cleared its natal molecular material, also preferen- 
tially located in regions of high extinction. 

In section [2] we describe the observations and compare 
several methods of selecting candidate YSOs. We will 
use the term YSO for any evolutionary state - heavily 
embedded, envelope dominated, disk dominated, or rem- 
nant disk (commonly called Class 0,1,11, and III for solar- 
mass objects). Sections |3] and m respectively describe the 
nature of the population and of individual YSOs in the 
region. The conclusions are summarized in fjs] 

2. OBSERVATIONS AND ANALYSIS 

M16 was observed as part of the Galactic Legacy In- 
frared Mid-Plane Survey Extraordinaire (GLIM PSE), a 
Spitzer Legacy program ([Benjamin et al.ll2003D. The re- 
gion was imag ed with IRAC (" Fazio et al.ll2004f ) at 3.6, 
4.5, 5.8 and 8.0 ^m with two 1.2s exposures at each lo- 
cation. The point source sensitivities are approximately 
0.5, 0.5, 1.0, and 4.0 mJy in the four filters, somewhat 
brighter in regions of intense diffuse emission. Sources 
were extracted from the images using the GLIMPSE 
pipeline, built around a mo dified version of the PSF- 
fitting program DAOPHOT (jStetsonlll987l Babler et al 
in preparation). Source detection and extraction is done 
in an iterative fashion. Initial sources were detected at a 
3-sigma level above the local background. These sources 
were subtracted from the image and then a second round 
of detection was run at a 5-sigma level to reveal sources 
that may have been missed in the wings of the PSF 
of bright sources. After the second round of sources 
were extracted a third iteration of photometry was done 
to improve the fluxes. In regions of complex diffuse 
background, the source extractor can sometimes over- 
estimate or under-estimate the point source flux (due 
to the difficulty of determining the local background). 



To remedy this situation, aperture photometry was done 
at the positions of extracted sources on the residual im- 
ages (the source subtracted images) to check for excessive 
over- or under-subtraction of the source. If the initial ex- 
traction was indeed over- or-under subtracted, then the 
aperture photometry was used to adjust the initial ex- 
traction value. The pipeline has been rigorously tested 
against other extraction software and using a network 
of calibrator stars with flux densities predicted by spec- 
tral type and multiwavelength photometry (M. Cohen, 
private communication). Our simulations indicate that 
the flux densities of truly pointlike sources are in gen- 
eral accurate to ~10% even in regions of high diffuse 
background. Zero point flux densities (Vega magnitude 
system) used in the GLIMPSE vl.O processing and this 
work are 280.9,179.7,115.0, and 64.13Jy for the IRAC 
bands, and we used a zero point magnitude of 7.16Jy 
for our MIPS 24jnm p hotometry (The IRAC calibration 
of iReach et al.ll2005L is slightly different, and users of 
version 2 or higher GLIMPSE products should refer to 
the accompanying documentation). Sources detected at 
several wavelengths are cross-correlated (bandmerged) 
using the SSC banderger software, which is based on 
the 2MASS bandmerger. Potential matches are consid- 
ered within 4" , and in confused regions, the association 
between all bands which minimizes the positional chi- 
squared is chosen. 

For this study, some of the sources of particular inter- 
est (e.g. in the tips of the pillars) were re-photometered 
by hand for this study, using irregular apertures and 
careful sky selection in regions of irregular diffuse emis- 
sion. The results of the manual photometry were for the 
most part consistent with the pipeline extractions, but 
in a few cases of slightly extended sources, the manual 
fluxes were substituted. It should be noted that photom- 
etry is extremely difficult in very crowded regions with 
structured diffuse emission such as M16, and although 
the GLIMPSE pipeline performs as well or better than 
any other extractor we know of for IRAC data, one still 
needs to be very careful about cross-correlating sources 
between filters and resolving close pairs of sources at one 
wavelength but not another. As we discuss below, exam- 
ination of each source's spectral energy distribution, and 
comparing it to model SEDs, is actually a very good way 
to weed out bad matches and questionably photometry. 
Sou rces were bandmerged with the 2MASS all-sky cata- 
log (jSkrutskie et al.ll2006( ) using the GLIMPSE pipeline. 

The GLIMPSE pipeline extracted over 42000 sources 
in the region (within 0.45 degree of ^=17.0, 6=0.8, near 
NGC 6611), and we wish to classify all of them in a con- 
sistent manner, as potential YSOs, or reddened main se- 
quence or post-niain-sequence stars. We considered and 
compared four methods: 1) NIR-excess sources in a JHK 
color-color plot 2) sources with a 2-24 //m spectral index 
greater than -2.0 3) sources in the IRAC color ranges of 
lAllen et all (l200l 4) sources better fit by a protostellar 
than stellar model SED. 

Figure [T] shows a NIR color-color diagram dissected 
into three regions in the usual manner. The two dividing 
lines are parallel to the r eddening vector. We use E(J- 
H)/E(H-K)=1.75 as in llndebetouw et al.l (|2005f ) and 
iRieke^fc Lebofsky (1985^, although a range between 1.6 
and 2.0 has been obser ved in different molecular clouds 
and the diffuse ISM (jCardelli et al.l Il989l : iFitzpatrickl 
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Fig. 1. — NIR color-color diagram with two lines parallel to the 
reddening vector. These separate the loci of dwarfs and giants, 
T-Tauri candidates, and protostar candidates (the least evolved 
young stellar objects, Class O/I for solar-mass YSOs), from top 
to bottom on the plot (see text and Meyer et al. 1997, for dis- 
cussion). Candidate YSOs selected by SED fitting are marked 
in orange (here and throughout, we will use the term "YSO" 
to includ e a ll evolutionary stages, not just the youngest). As 
noted in [ Ro bitaille et al. (2006) and shown observationally (e.g. 
IMaercker et a l. 2006; Strom et al. 1989 ). only the small fraction of 
YSOs with unobscured very hot dust have significant JHK excess. 

119991 ). The line originating at (0,0) separates reddened 
main-sequence and giant stars from T-Tauri candidates. 
The other hue separates T-Tauri candidates from less 
evolved YSO candidates (a dotted line marks the classi- 
cal T-Tauri locus, or experimentally determined color lo- 
cus of unredenned T-Tauri stars - disk-dominated YSOs 
are expected to lie near this line or along a redenning 
vector originating on the line iKenvon fc Hartmannlll99"5t 
iMever et aLlll997D . We define the "NIR excess" as the 
distance down and to the right of the first line (reddened 
MS and giant stars have excess < 0, YSO candidates 
have excess > 0.3): 

{H ~ K)€ + 0.59 ~ {J ~ H) 
excess = ^ — , — ^ ( 1 ) 



e^E{J - H)/E{H - K) 



(2) 



Most other classification schemes for young stellar ob- 
jects use mid-infrared colors, or the mid-infrared spectral 
index. Classically, the 2-24 fim spectral index (power- 
law exponent fitting AFa oc A") has been used to classify 
YSOs, and to define the "Class I-II-III" system (Lada, 
1987). In the Spitzer era, many star formation regions 



have been mapped with arcsecond spatial resolution us- 
ing the 3.6-8.0 /im IRAC camera (|Fazio et al.ll2004l ). so 
it is useful to consider methods of selecting YSOs based 
in IRAC photometry alone. However, Whitney et al. 
(2003a, b, 2004b) showed that the narrower wavelength 
range (3.6-8.0 /im v. 2-24 ^m ) is more susceptible to 
2-D effects (inclination), and stellar temperature. The 
8.0 fiui band is also affected by the broad 10 fim silicate 
feature that could be either in emission or absorption. 
More recently, we have computed a large grid of YSO 
models and found that while the 3.6-8.0 /.tm spectral in- 
dex does correlate roughly with evolutionary state (mea- 
sured by circumstellar dus t mass and accretioii rate), 
there is significant spread (jRobitaille et al.l [20061 Fig- 
ure 11). 

Classification based on IRAC colors uses more infor- 
mation about the SED than a simple spectral index over 
the same wavelength range, and ma y be more robust . 
Figure [2] shows the regions defined bvlAllen et all (l200j, 
hereafter A04) (see also iMegeath et all \2004) in IRAC 
color space. Based on t he solar-mass protostellar mod- 
els of lCalvet et alj (|1994[ ). those authors suggest that the 
sources in the central box are likely to be Class II YSOs, 
and redder sources Class I or younger. Our models do 
not show as neat a separation in color-color space as A04 
due to th e larger range of para meter space explored in 
our grid teobitaille et al.ll2006l ). The deeply embedded 
sources show color variations due to inclination effects 
from scattering in the bipolar cavities, and sources of 
all evolutionary stages show color variations due to inner 
hole size and stellar temperature (Whitney et al. 2004b). 
In Figure [21 we compare A04 IRAC color-selected proto- 
stellar candidates with the spectral index and the NIR 
excess (for sources with J,H, and Kg detections). The 
IRAC-color classification is very consistent with having 
large 2-24 fim spectral index, and somewhat consistent 
with having a NIR excess, although the latter is not ex- 
pected for all young sources. 

In this study we use a different type o f classification of 
YSO candidates - model SED fitting (jRobitaille et al.l 
|2007| ). This is accomplished using a program that 
fits all available fluxe s to a large grid of model SEDs 
(jRobitaille et al.l [20071 here we use fluxes between 1.25 
and 24 /xm, but any between the UV and millimeter are 
usable). Our model grid consists of 2x10^ spectral en- 
ergy distributions emitted by 2x 10^ young stellar objects 
(YSOs; 10 inclination angles each). The model grid is de- 
signed to span both theoretically and observationally rea- 
sonable ranges of mass, temperature, evolutionary state, 
disk mass, disk accretion rate, envelope mass, and en- 
velope accretion rate. The colors of the model grid also 
span the full range of observed colors of stellar and young 
stellar objects. In addition to the 200000 YSO models, 
our grid also includes 78 53 stellar model atmospheres 
(jBrott fc Hauschildtll2005h and 75 AGB star templates 
derived from ISO spectroscopy and spanning the range 
of evolutionary state and mass loss rate. The fitter re- 
turns a best-fit for each source as well as the distribution 
of templates and models that fit with smaller than a user- 
selected value of ■ We first select all sources that can 
be fit well (x^ per datapoint, or xt < 2-5) by an extin- 
guished stellar atmosphere, and consider these consistent 
with stars. Of the remaining sources, if one can be fit 
with xt < 2.5 by any YSO models, then we consider it a 
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Fig. 2. — The first panel shows the IRAC color-selection areas 
defined in A04,as applied to sources in the M16 region. The central 
box encloses candidate Class II sources (orange), and the redder 
sources candidate Class I. The grey points are assumed to be stars 
with no circumstellar dust. An extinction vect or for Ay =30 using 
the measurements of Indebetou w et al.l l|2005l ) is shown. The sec- 
ond plot shows the A04 classification compared to the NIR color 
excess and MIR spectral index. Sources at the right edge of the 
plot do not have J-band detections so no NIR excess could be cal- 
culated. Sources at the bottom do not have measurements at 2 or 
24/im, so that particular spectral index cannot be calculated. The 
right-hand axis shows the spectral index that a r eddened Rayleigh- 
Jeans photosphere would have at a given Aj^ IHillenbrand et al.l 
[1993. measured Ay = 3.2 towards NGC6611). 

candidate YSO. If it cannot be fit well by either YSO or 
stellar atmosphere, then it is likely poor photometry or 
bandmerging. This was confirmed by manual examina- 
tion of the SEDs, and the process is robust to changes of 
the threshold goodness-of-fit between 2.0< xt <4.0 (the 
number of sources in each category only changes by a 
few percent). For this analysis we consider YSO models 
with distances betwee n 1.5 and 2.5 kpc - as discussed in 
iRobitaille et al.l (|2007f ). if the distance to the sources can 
be constrained, the models can be better constrained by 
the data. 

Analyzing as much of the SED as possible has sev- 



eral advantages. In crowded regions with bright com- 
plex diffuse emission such as M16, examination of the 
SED provides an extra check to remove bad data points. 
For example, we found several sources that were resolved 
from a close companion at H band and by IRAC, but not 
in the 2MASS data, resulting in anomalously high 
Kg fluxes (which would be interpreted as a NIR excess). 
SED analysis is robust to instrumental artifacts and bad 
data in a single filter - color-color analysis requires good 
photometry in four filters (the four IRAC filters for ex- 
ample) , whereas we can admit measurements with larger 
uncertainties in some filters, but are protected against in- 
troducing many false YSOs because if the resulting SED 
is inconsistent with any stellar or protostellar model, it 
can be eliminated from the analysis. 

Perhaps most interesting is the ability to understand 
the constraints afforded by the data in as objective a way 
as possible. If the photometry of a given source is con- 
sistent with a wide range of circumstellar dust models, 
but inconsistent with any stellar atmospheres, the source 
can still be classified as a probable pre-main-sequence ob- 
ject, but at the same time it is known that the data do 
not allow the evolutionary state (disk-only or accreting 
envelope) to be determined. This is philosophically sim- 
ilar to choosing all sources with spectral index greater 
than some threshold, or all sources with some threshold 
infrared excess. In this manner, we use the model fit- 
ter to determine upper limits to the evolutionary state of 
sources in M16: Namely, if a source is not fit by any stel- 
lar atmosphere with xt <2.5, but is fit at least that well 
by some YSO models, the source is considered a YSO 
candidate. If the source is fit by models with accreting 
envelopes, but cannot be well-fit by any disk-only model 
(no accreting envelope), then we consider the source a 
candidate to be quite young (perhaps Class I for a solar- 
mass object). Finally, if the source can be fit by mod- 
els with circumstellar disks more massive than 10~^xM* 
(where M^, is the central source mass in the model), but 
is not fit by any models with less massive disks, then 
we consider the source to be a candidate accretion disk 
source (similar to Class II for a solar- mass object). 

Figure [3] compares the classification suggested by our 
model fitter with the NIR excess (for sources which JHK 
detection allows calculation) and MIR spectral index. 
Red points have a high probability of being relatively 
unevolved with massive envelopes, orange points likely 
have still-accreting circumstellar disks, and green points 
are YSO candidates for which the existing data cannot 
place strong constraints on the circumstellar dust distri- 
bution. We see that YSO candidates generally fall in 
the regions defined by A04. The constraints suggested 
by the model fitter are less well-defined than the box 
drawn by A04, as expected given a range of inclination, 
stellar masses, and inner disk holes. One should also re- 
member that the colors plotted are upper limits to the 
evolutionary state, and the source may have more mas- 
sive envelopes or disks, but the existing data do not allow 
that to be unambiguously determined. Foreground ex- 
tinct ion_afsojnovBS_souree^ in this color space 
(e.g. Ilndebetouw et al.|l2005[ ). which explains some or- 
ange points lying above the box. Some of the orange 
points to the right of the box (large [5.8]-[8.0], moderate 
[3. 6] -[4. 5]) have inner holes in their dust envelopes. In 
the case of sources for which the data can constrain the 
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circumstellar dust distribution, there is a correlation be- 
tween less evolved objects (red points) and higher MIR 
spectral indices. 
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Fig. 3. — Comparison of classification using our model fitter 
with the color regions of A04 (first panel), and MIR, spectral index 
and NIR excess (second panel) . Axes and annotations are identical 
to Figure [2] An extinction vecto r for Ay=30 using the measure- 
ments of Indebetouw et al] II2005I ) is shown on the first panel. Red 
points have a high probability of being relatively unevolved with 
massive envelopes, orange points likely have still-accreting circum- 
stellar disks, and green points are YSO candidates for which the 
existing data cannot place strong constraints on the circumstellar 
dust distribution, but which are better-fit by YSO models than 
stellar atmospheres (see text for details). The sources may have 
more massive envelopes or disks than shown, but the existing data 
do not allow that to be unambiguously determined. 

Several characteristics of the fitter classification com- 
pared to a classification based only on IRAC color are 
particularly important. Comparing the second panels of 
Figures and [31 we see that the multi-band fitter finds 
many more young sources that did not have JHK detec- 
tions (the dots on the right-hand side of Figure [3]). Our 
fitter can characterize a source with measured fluxes in 
as few as four bands between J and [24] , whereas a NIR 



excess characterization clearly requires JHK detections 
and an IRAC-color characterization requires detections 
in all four IRAC bands. The latter can be difficult in 
regions of high diffuse emission, like a star formation re- 
gion. (A similar effect to fitting an SED to all bands 
can probably be achieved by classifying sources in mul- 
tiple color-color planes; Megeath et al. 2005, Allen et 
al. 2005). Careful examination of the entire SED al- 
lows the removal of confused sources or those with poor 
photometry, which is less obvious in a color-color (NIR 
or IRAC) plot. Finally, comparison of the SED with as 
many different YSO models as possible provides a rela- 
tively unbiased way to assess the degree to which the data 
can constrain the modeled nature of the source, beyond 
assigning a most likely model based on colors. 

3. THE DISTRIBUTION OF YSOS IN THE EAGLE NEBULA 

Figure [4] shows the Eagle Nebula in the [3.6], [4.5], and 
[8.0] bands. The optical open cluster is the concentration 
of blue stars to the NW of the HST pillars (marked in 
Fig [5]). Protostellar candidates selected using our fitter 
are shown on a [8.0] image in Figure [5l Many of the pro- 
tostellar candidates are located in the NGC 6611 cluster, 
which is known to be very young. In fact, nearly all of 
the sources in the NGC 6611 region are classified as YSO 
candidates: 46 of the 47 sources with high signal-to-noise 
detections in at least 4 bands between 1 and 10 fim and 
located in a 4' diameter region around NGC 6611 were 
classified as YSO candidates. This is consistent with a 
very high L- band (3.5 iim) exces s fraction (85±5%) de- 
termined bv lOliveira et al.| (l2005f) and the large nu mber 
of emission-line stars seen bv lde Winter ir^|l99l . im- 
plying that circumstellar disks are fairly robust, even in 
a cluster environment. 

Our analysis reveals a new cluster to the NE (at 
18hl9ml0s, -13d36m, approximately 7' in diameter). 
Located in a dark region indicative of a dense cloud, 
the new cluster is likely younger than NGC 6611 - as 
seen in Figure [5l most of the new cluster's sources can 
be constrained to have accreting disks (yellow) or en- 
velopes (red), compared to NGC 6611 which has many 
sources consistent with remnant disks (green). There 
is also a concentration at the head of a large pillar at 
18hl9m05s -13d45m20s ("Pillar 5"; we analyze this re- 
gion in more detail below, H4.3p . Finally, there is a low 
density of scattered protostellar candidates across the re- 
gion shown. This apparent interspersed ongoing star for- 
mati on (also n o ted in the giant HII reg ion RCW49 by 
Whit ney et all ()2004aD and in Orion by iMegeath et~al] 
could possibly be a superposition of background 
objects, but at least some of them are likely YSOs, 
indicating widespread low-level activity throughout the 
cloud. 

3.1. Spatial Relationship of YSOs to the Molecular 

Cloud 

A byproduct of our fitting program is the best-fit ex- 
tinction to each source, including the large number of ex- 
tinguished main-sequence and giant stars along the line- 
of-sight. These values can be binned to produce a map 
of the average extinction, which is shown as contours 
in Figure [HI We expect that the molecular cloud associ- 
ated with star formation in M16 dominates the extinction 
measured with this procedure, and indeed the features 



6 



R. IndebetoTiw et al. 




Fig. 4.— . The Eagle Nebula as seen by IRAC (RGB 
NW corner. 



[8.0], [4.5], [3. 6]). IRAC [8.0] has been supplemented with MSX A (8 /im) in the 



match other indicators of dense materiaL NGC 6611 is 
located in a low-extinction region, likely evacuated by the 
cluster itself. High-extinction regions correspond reason- 
ably well with the boundaries of the bubble delineated 
by PAH emission in the the IRAC bands and scattered 
light in the NIR. There is a concentration of dense gas 
to the NE of NGC 6611, and more star formation on 
the other (north) side of the cloud, in regions of some- 
what lower extinction. We compared the spatial distri- 
butions of YSOs as a function of evolutionary state and 
mass, and found no obvious variations with either pa- 
rameter other than the likely younger NE cluster. Since 
our maximum extinction through the cloud is Ay —20, 
which is 1 magnitude of extinction at IRAC wavelengths 
(jlndebetouw et al1l2005D . we can still detect YSOs down 
to masses of about O.SIMq, and we do not expect to be 
missing any significant star formation in the region. 

The combination of mid-IR imaging at different wave- 
lengths allows one to directly compare the protostellar 
population to the properties of the molecular cloud in 



which they are forming. First, we use IRAS 60 /im and 
100 fim images to determine the temperature of the dust: 



Tdust — 



hc/k(l/\2-l/\l) 



(3 + (3) ln(Ai/A2) + ln(F„(Ai)/F,(A2)) ' 
where (3 is the emissivity index, 



2hc 



A 



;^3(-g/tc/(AfcTd„,t) _ 

(see e.g. lUfc Chenlfl99l iSchnee et al.l[2005l) . Figure[7] 

shows contours of the dust temperature relative to the lo- 
cations of YSOs and PAH emission. Clearly, the warmest 
dust is heated by NGC 6611, as one would expect. It is 
also evident by the NE and SW extensions of the con- 
tour shapes centered on the "pillars of creation" that 
the cavity presumably blown by that young open cluster 
is warmer than the surrounding denser molecular cloud. 
There is also some evidence of heating to the west of 
NGC 6611, perhaps by embedded sources like the mas- 
sive YSO that we find there 
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Fig. 5. — IRAC [8.0] image of M16 with stars of various inferred evolutionary states marked in different colors. Sources that can be 
constrained to have relatively massive accretion disks are yellow, and those that in addition likely have massive circumstellar envelopes are 
red. YSO candidates for which the data do not strongly constrain the mass of circumstellar dust (merely that they have some, and are 
inconsistent with stellar atmospsheres) are green. NGC 6611 and a new cluster of YSO candidates to the north are marked with large blue 
circles. The edge of the region for which coverage exists in all IRAC bands and with MIPS at 24 /im is marked as a thin black line. 



The optical depth of a molecular cloud can also be cal- 
culated from the dust emission, and this compares well 
with our map calculated using embedded stars. Given 
a temperature, the optical depth of the emitting dust 
is simply Ti, = / By{Tdiist), where B^{T) is the Plank 
function. We use the 100 /im image to calculate t, but 
the results using the 60 /.tm are very similar. We used 
(3 =2 for these calculations; choosing (3 =1 increases the 
temperatures by 15% and decreases the optical depths 
by a factor of 2, but does not affect the relative optical 
depths, which features we are most interested in exam- 
ining here. Figure [8] compares the optical depth calcu- 
lated by the fitter using point sources photometered by 
GLIMPSE to the optical depth calculated from cool dust 
emission - most features agree well, including the cavity 
and dense ridge. The new NE cluster region shows high 
extinction in the point sources, but low optical depth in 
the dust. 



3.2. Clustering 

Stars are generally thought to form in cluster envi- 
ronments, so greater clustering in a population should 
be a sign of youth. The clustering properties of young 
stellar objects should also reflect any preferred spatial 
scales associated with the formation mechanism, if those 
exist. In Figure [9] we show the two-point correlation 
function of our YSO candidates in the M16 region. We 
count the number of pairs of YSOs as a function of an- 
gular separation 0, binned logarithmically - i.e. N{6) 
is the number of pairs with separation between 9 and 
9 + d\og9. This distribution is compared to the dis- 
tribution of separations of a random field of stars with 
the same average stellar density Nr{9); the correlation 
function is N/Nr — 1. The maximum value of the corre- 
lation function is comparable to peak values of 0.5, 1.5, 
6, and 15 in subcl usters in the Rosette Molecular Cloud 
(|Li fc Smithll2005l ). reflecting the strong degree of clus- 
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Fig. 6. — [8.0] fim image of M16 with contours of average ex- 
tinction and the location of protostellar candidates (color as FiglSj. 
Contour levels are Av=10.5, 12.5, 14.5, 16.5, 18.5, 20.5. NGC 6611 
is located in a low-extinction region, likely evacuated by the cluster 
itself. There is a concentration of dense gas to the NE of NGC 6611, 
and more star formation on the other (N) side of the cloud, in a 
region of somewhat lower extinction. The dark area to the NW 
is missing IRAC data, so the blue contours end artificially due to 
lack of point source photometry 




Fig. 7. — Contours of dust temperature calculated from IRAS 
60 /jm and 100 /im images, overlaid on the IRAC [8.0] image of 
M16. Contours levels are 32, 33.3, 34.7, 36, 37.3, 38.7, and 40K. 

tering. The correlation function is only roughly consis- 
tent with a power-law (Fig[9l second panel), which would 
mean there is no preferred spatial scale for star forma- 
tion in the region. The slope of the correlation function 
in logarit hmic space (-0-8 ) is also comparable to other 
clusters (Li fc Smithll2n05l : iScalo fc ChaDDelll[l999l ). also 
indicating fairly strong hierarchical clustering. Typical 
optical depths in region are Ay=15-20, corresponding 
to N(H2)^4xlO^*' in the molecular material. The cloud 
appears filamentary, with thickness about 4pc, so if the 
depth along the line of sight is similar, the average vol- 
ume density in the cloud is n(iJ2)'~^150cm~'^, with a ther- 




FlG. 8. — Contours of extinction calculated by fitting 42000 stellar 
point source SEDs, overlaid on a gray scale image of the optical 
depth of emitting dust calculated from IRAS 60 fim and 100 fim 
images. Candidate YSOs identified by our fitter are marked as 
circles colored as Fig[5] The dark area to the NW is missing IRAC 
data, so the blue contours end artificially due to lack of point source 
photometry. 

mal Jean's length of 2pc (Rj ~4.8pc(T3/n)° '5 if T is ex- 
pressed in K and n in cm~'^). At a distance of 2kpc, 
this corresponds to a size scale of 10"^ '^; there may be 
some signature of a turnover in the correlation function 
at about 10"^ '' degrees, but the result is not definitive. 
There are also structures in the region of about that size 
(for example the bubble around NGC 6611, and star- 
forming pillar 5, f ijlip and the correlation function could 
be reflecting those structures. 

The distribution of nearest neighbor distances of young 
stellar objects can potentially reflect a fragmentation 
scale. We show that distribution in Figure [TOl com- 
pared to a random distribution. The YSOs show an ex- 
cess at separations ^^10", or O.lpc at a distance of 2kpc. 
This is consistent with the fragmentation scale in the 
hot dense gas such as in the heads of the pillars and 
molecular cores near HH2 16, which are warm and dense 
(T-40K, n^lO'^cm- MPoundlfTOOSHAndersen et al.ll2004l . 
yielding Rj ~0.1pc). The signature of the parent molec- 
ular cloud's Jeans length has b een seen in the sepa ra- 
tion of YSOs before, notably bv lTeixeira et al.l (|2006f ) in 
NGC2264. The YSOs with nearest neighbor distances 
in the range 10~^'^ to 10~^'^ degrees are located in the 
densest filamentary parts of the molecular cloud to the 
north (seen as dark lanes in the 8 /um image) , associated 
with the HST pillars and pillar 5, and also in NGC6611. 
This association with dense molecular material suggests 
that the separation and Jeans length may not merely be a 
coincidence. It is less clear why many stars in NGC6611 
have that nearest neighbor separation - perhaps condi- 
tions in the cloud at the time of formation of NGC6611 
were similar. 

4. SPECIFIC OBJECTS AND REGIONS 

4.1. The HST "Pillars of Creation" 

Using ISOCAJVI, iPilbratt et "all ()1998[ ) stated that the 
lack of mid-infrared emission along the length of the 
pillars indicated that there was not enough dense star- 
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Fig. 9. — The two-point correlation function of YSO candidates 
in M16. The correlation function is plotted on a linear scale in the 
first panel, and on a logarithmic scale in the second panel. The 
dashed line is the correlation function of only the less evolved YSO 
candidates (best-fitting models have dust envelopes v^ith nonzero 
accretion rates). The dotted line on the first panel is the correla- 
tion function for all sources in the region, which are uncorrelated. 
The dotted line in the second panel corresponds to a power-law 
correlation function with index -0.8 



Fig. 10. — The distribution of nearest neighbor distances of 
YSO candidates in M16 (solid line) compared to the distribution 
for a random star field (dotted line). The YSOs show an excess at 
separations ~10", or O.lpc at a distance of 2kpc. 




■i 



Fig. 11. — Individual sources in the HST pillars region, shown 
on IRAC [8.0] image. See text for discussion of their nature, and 
Tables [T] and |2] for coordinates and physical parameters. 



forming material present to form new stars. The resolu- 
tion and sensitivity of the ISOCAM images was not suf- 
ficient to rule out very embedded or low-mass stars (the 
sensitivity limit of 20m Jy at 6.7 /im in bright regions cor- 
responds to a 4M0 YSO, Av=30, 10^ <age<10V)- The 
higher resolution and sensitivity of the IRAC images al- 
lows us to confirm Pilbratt's statement with much more 
stringent upper limits - we detect no new mid-infrared 
point sources down to a limit of '--^ImJy in the [3.6] band, 
and ~5mJy in the [8.0] band (limited by the strong dif- 
fuse background emission), corresponding to 0.4Mq. 

There are a handful of YSOs in the pillars, al- 
ready identified in t he NIR (e.g. ISugitani et"aLl 120021 : 
I Andersen et~ani2004t [Thompson et al.ll2002f ). Figure HH 
shows the locations and names of notable objects. Our 
new mid-IR flux densities for many of these allows a more 
precise determination of their likely physical parameters 



than was possible from NIR data alone. We use the 
grid of YSO models described in fJ^J and demonstrate 
how to make a careful assessment of all the models that 
are consistent with the data. Particularly important are 
any age constraints that could strengthen or weaken ar- 
guments that the YSOs in the trunks are triggered by 
the action of NGC 6611. The evaporation timescale for 
the trunks is 2xl0^yr, but internal velocity gradients 
sugges t that they w ill disperse or be torn apart in only 
lO^r (|Poundlll998l ). The trunks appear to be in rough 
pressure equilibrium with the sur rounding H II region 
(P/k~8xlOW-3K, lPoundlfl998h . so there is no well- 
defined "crushing" timescale. 

The most prominent young objects in the pillar re- 
gion are PI and Tl at the tips of the first and 
second most northern pillars (Hi and 112), respec- 
tively. At our modest 1.5" spatial resolution, the 
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PI sour ce is the same as that identified in the 
NIR, bv iSugitani et all (12001. and as M16ES-1 by 
iThonipson et aLl 1 2002[ ) . iFukuda. Hanawa. k. Sugitanil 
(I2002D detect 2.7mni sources further behind PI in the 
pillar and suggest that these are distinct YSOs, but our 
data does not support that suggestion. Figure fT2l shows 
the MIR SED and protostellar models consistent with 
the data. The 2-24 /im spectral index of t he source is 
a=l±0.2, Class I in the usual nomenclature ()Lada 1119871 
Class I: 0< a <3, Class II: -2< a <0, Class III: -3< a <- 
2). Source Pi's bolometric luminosity (estimated by 
YivF^) is 44±1OL0 (observed quantities are tabulated 
in Table jH). Using the pr e-main-sequence evolutionary 
tracks of lSiess et al.l ()2000() . this is consistent with a 2.5- 
4Mo 3x10^ old YSO (the tracks are quite close together 
so a small uncertainty in age or luminosity leads to a large 
uncertainty in mass). The YSO model fitter produces 
similar results, as tabulated in Table [H Models with 
masses between 1.6 and 6.4M0 fit the data with xt <3.7, 
and the best-fitting model has a mass of 4.5M0(for this 
analysis we consider all models which fit to <3.7 be- 
cause that represents a probability of 1/1000, and we 
wish to be conservative in our constraints on physical 
parameters) . Those models fit with a foreground (line-of- 
sight) extinction Ay ~4. They have envelope accretion 
rates M between 1 x 10~^ and 7xl0~'^M(T)/yr, also consis- 
tent with a Class I source. [Thompson et al.l ()2002l ) derive 
a stellar luminosity of 2OOL0, and comment that the lack 
of Paa emission suggests a group of several lower-mass 
objects. Our fitted n iass is that of a very late B star 
and low ionizing flux. iMcCaughrean fc AndersenI ()2002f ) 
find a somewhat higher mass (IOM0, a Bl star). Those 
authors also derive an extinction of Av=21 , if all the J 
and H hand light is photospheric. The extinction from 
our model fit of Ay —4 listed in the Table is only the 
foreground extinction to the YSO; the total extinction 
including foreground and circumstellar d ust is Ay 0^40 

Ano ther red source, which we and ISugitani et al.l 
(j2002f) call P2, is located (at least in projection) below 
the tip of Pillar 1 with respe ct to NGC 6611, and noted 
to ha ve high extinction by IMcCaughrean &: AndersenI 
(|2002|) . who call it E23 in correspondence with the H96 
EGG of that designation. This source has a spectral in- 
dex of 0.3±0.2 and a bolometric luminosity of 5.3±2.8L0 
- the inability to determine a reliable 24 /im flux makes 
it more difficult to constrain this object's luminosity and 
evolutionary state. Using the PMS tracks as for PI and 
assuming an age of a few 10^ years suggests that this 
source is about IM0. Using only a 24 /xm upper limit 
does not allow the model fitter to constrain the param- 
eters very well either, with accretion rates between zero 
and 10"'^ and masses between a few tenths and a few 
solar masses fitting the data with xl < 3.7. The mid- 
infrared excess and positive spectral index suggest that 
accretion is still continuing in this object (and no model 
without any accretion disk is consistent with the data). 

The source at the tip of the central pillar (Tl) was clas- 
sified as a somewhat olde r source, a candidat e T-Tauri 
star, based on NIR excess (ISugitani et al.ll200l . We find 
a spectral index of -0.3±0.1 and luminosity of 22±8L0. 
This lower spectral index falls in the Class II (T-Tauri) 
range, and a several 10^ year old object of this lumi- 
nosity falls in the 2-3M0 mass range of the ISiess et al] 
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Fig. 12. — Source PI. The first panel shows the data and the 
model SEDs that fit better than xi < The second panel 

shows the envelope accretion rate M and disk mass M^ji^^; for 
those models. The pl ot range is that of the entire model grid (see 
IRobitaille et al.1 120061 . for details). The bars along the axes are 
for models with zero M and M^j^j., i.e. if there were well-fitting 
models with zero disk mass they would be in the bar along the x 
axis (this will be true for some sources below). Nearly all of the 
models that fit the data for this source well have high accretion 
rates and disk masses, consistent with the large spectral index and 
a relatively young evolutionary state. 



(|2000f ) PMS tracks. The models consistent with the data 
have similar masses (2.4-5.3M0), and relatively low ac- 
cretion rates (many disk-only models, without an accret- 
ing envelope, fit the d ata well as seen in Figure [T4| . 
iThompson et al.l (|2002| ) derive a rather lower luminosity 
but had to extrapolate from low -resolution MIR data. 
IMcCaughrean fc AndersenI (|2002f) determine a compara- 
ble mass to ours, 2-5M0. They determine an extinc- 
tion "down to the photosphere" of Ay = 15; we deter- 
mine a foreground extinction of ^--^5, and the models that 
fit with xt < 3.7 have a range of circumstellar extinction 
from <1 to ~20. As seen in Figure [TH the SED would 
be better defined, and the circumstellar dust distribu- 
tion better constrained, with optical measurements or a 
more precise 24 /zm flux density (which would require 
higher spatial resolution to better separate the source 
from nearby diffuse emission from the pillars). If indeed 
this object is fairly evolved (several million years old), 
then it seems unlikely that it was triggered by the ac- 
tions of NGC 6611, unless there really are 6Myr evolved 
stars in that open cluster that were able to make their 
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Fig. 13. — SED and model fits to the source P2. Plots and 
annotation as Fig |12l 

effects felt very early. It is also possible that the pho- 
toablating effects of NGC 6611 have caused premature 
dispersal of the larger envelope around this source. 

The northernmost pillar (Hi ; Fig [TT|) is actua lly two 
head-tail structures, or two molecular clouds (Poundl 
Il998f ). YSOs found half-way down the pillar obey 
the same paradigm of "sources forming at the tips 
of pillars" as the sources at the three tips. After 
YSOl(Pl) and YS02(T1), two of the objects with 
highest estimated mass and extinction in the pillars 
are located in this region: H96 EGGs E31 and E42 
(jMcCaughrean fc Andersenir2002 ). We wi ll refer to these 
as P3 and T2, respectively, following ISugitani et al.l 
(|2002f) . P3 (E31) is one of the brightest objects in the 
pillars by 8 /im, much brighter than its optically-bright 
neighbor to the SW. It has a very positive spectral index 
2.2±0.2 and luminosity of 15±5L0, which corresponds 
to about a 1.5Mq 2x10^ year old PMS star. Like PI, 
only models with relatively massive, accreting envelopes 
fit this source well; it must also be quite young. In both 
sources, we see that when an object is modeled as hav- 
ing a massive envelope, it is more difficult to constrain 
the mass of the accretion disk inside that envelope from 
infrared data alone. 

Source T2 is similar in nature to Tl: its spectral in- 
dex 0.2±0.2 is consistent with zero, and the majority of 
well- fitting models have low accretion rate (see Fig [TB)) . 
Its fairly low luminosity 6±1.5Lq corresponds to l.SM© 
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F ig. 14. — Model fits to the source Tl. Plots and annotation as 
Fig [12] 

at 10^ years, and although the best-fitting models have 
somewhat higher masses, the lack of optical and 24 /im 
fluxes make it difficult to constrain the parameters of 
this source beyond stating that the probability of a very 
massive envelope is low. Again, if this object is truly so 
evolved, then the case for triggering by l-2Myr old stars 
in NGC 6611 is weakened. 

Two stars are located in the concave tip of Ila, ei- 
ther slightly more evo lved objects, or a chan ce line-of- 
sight alignment (H96. [Thompson et al.ll2002f ). We find 
that these objects are best-fit with extinguished main- 
sequence stars, not YSOs, supporting the theory that 
they are not very young objects. Other previously noted 
objects in the region, such as sources "I" (to th e north 
of Pillar 1) and "367B" of iPilbratt et all (|1998D (to the 
east of P3) are also best fit by extinguished photospheres, 
and we do not consider these strong candidates for young 
stellar objects. 

4.2. HH216 

IIII216 is an optically visible Herbig-Haro object lo- 
cated just sou th of the base o f the three HST pillars 
(see Fig [TT])- lAndersen et alJ ()2004f ) propose that the 
driving source of HH216 is one of two very red sources 
located in a structure to the south of the HST pillars, 
which they call pillar #4. The two sources are irregu- 
larly shaped in those author's NIR data, and are clearly 
detected in the IRAC data, although unresolved with our 
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F ig. 15. — Model fits to the source P3. Plots and annotation as 
Fig [12] 

poorer resolution. In fact, the confused diffuse emission 
in the region makes it difficult to extract a reliable flux at 
[8.0], leaving the SEDs less constrained than some of the 
other objects which we have discussed. The shape of the 
SED suggests quite high extinction {Ay of several tens), 
and the spectral indices of 1.6±0.6 and 1.1±0.5 for the 
northern and southern source respectively suggest youth 
and/or very high extinction. There is not clearly any 
associated 24 /im emission, ruling out extremely high lu- 
minosities. The Spitzer data are in all consistent with 
the solar-mass Class I sources that most commonly drive 
HH jets. Our data do not suggest one or the other as the 
more likely source for HH216. The source apparently 
located more in the tip of the pillar, at least in projec- 
tio n, which we call HH-N, has mas er emission detected 
bv lHealv. Hester, fc ClaussenI ()2004f ). a possible indicator 
of greater activity. 

4.3. Pillar 5 

There are several YSO candidates in a prominent pil- 
lar to the east of the HST pillars, near 18hl9m07s 
-13h45m24s. This pillar was called pillar #5 by 
iHealv. Hester, fc ClaussenI (|2004D , who discovered water 
maser emission in several spots. Figure [T7] shows the 
IRAC image of pillar 5. The masers are coincident with 
a bright red MIR source; in fact the source is marginally 
resolved in our images, with two components correspond- 
ing to two of the three different clusters of maser sources. 
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F ig. 16. — Model fits to the source T2. Plots and annotation as 
Fig [12] 

We call this source "A", and its SED is shown in Fig- 
ure [TH] This source's high spectral index 2.0±0.1, high 
luminosity 250±100Mq, and the fact that it is a bright 
24 /im source all make it a good candidate for a young 
intermediate mass YSO. At an age of l-2x lO^yr, that lu- 
minosity corresponds to ~6Mq . Well- fitting models have 
relatively high foreground extinctions Ay ~7 and quite 
massive envelopes of a few 1O~^M0. The other sources 
in the tips of the "fingers" of this pillar are also classified 
as YSO candidates, mostly with modest accretion rates. 
Another quite luminous young object is located near the 
base of the pillar, and is detected easily at 24 /im and 
even at 70 /im . Source "B" has a high spectral index 
1.7±0.1 and luminosity 2OO±5OL0. 

4.4. North Bay 

Another region highlighted by 

IHealv. Hester, fc Claussel] (|2004D for the presence 
of water maser emission is to the north of the HST 
pillars, around 18hl8m46s -13d44m30s. There is a 
bright submm source here, a faint green MIR source at 
the peak of the 850 fxin emission (not quite coincident 
with a water maser), and another offset MIR source 
which is exactly coincident with another water maser. 
The former, faint green source is only detected at [3.6] 
and [4.5]. The latter, on the ridge, is detected only 
longward of [3.6], and is only well- resolved from the 
diffuse emission at [4.5]. With so few extractable fluxes. 
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Fig. 17. — IRAQ [8.0] image of eastern pillar 5. Positions of water 
masers from lHealv. Hester, fc ClaussenI pOOS) are coincident with 
source "A". 
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Fig. 18. — M odel fits to the pillar 5 source A. Plots and annota- 
tion as Fig |12l 

it is difficult to say anything about the nature of these 
sources, except for their [3.6] and [4.5] flux densities: 
The faint green source has flux densities of l±0.5mJy 
and 10±3mJy, and the brighter source on the ridge 
13±5mJy and 20±3niJy in those two bands. 



Fig. 19. — IRAC image of North Bay region. SCUBA 850 /im 
contours are superimposed (data from the SCUBA archive, ob- 
servation date Feb 18, 2002, program M01BU2 0), and locations 
of wa ter masers marked as diamonds (jHealv. Hester, fc ClaussenI 
[200l) . 

4.5. a neAii Massive Young Stellar Object 

The most luminous YSO which we identify in the re- 
gion is labeled "MYSO" in Tables [1] and [21 and shown in 
Figure[20l This YSO is an IRAS source, 18152-1346, and 
it also appears to have associated water maser emission 
(jBraz fc E pchtein 1 983 ) , another sign of massive star for- 
mation. The shorter wavelengths (2MASS, IRAC [3.6], 
[4.5]) show a bright central point source surrounded by 
either a small cluster, not resolved into individual stars, 
or extended diffuse emission as might arise in an out- 
flow. The IRAC colors do not show a particular excess 
at [4.5], as sometimes occurs in sources believed to have 
outflows (and emission in that band from CO and II2 
lines). Precise photometry of this source is made diffi- 
cult by the extended nature at short wavelengths, and 
that it is saturated at [8.0] and [24], and marginally sat- 
urated at [5.8]. We used a combination of PSF wing 
fitting, aperture photometry, and comparison with the 
MSX [8] and [21] fluxes to derive the source's SED shown 
in Figure [211 The source's spectral index 1±0.3 and lu- 
minosity 1000±300Lq suggest that it is young, and may 
have a mass (from PMS tracks) of about 8Mq. Exami- 
nation and fitting of the SED shows likely PAH emission 
from the source, increasing the fluxes at [5.8] and [8.0]. 
PAH emission is expected from the surface of the accre- 
tion disk, or interior of an outflow cavity in a flattened 
envelope, for such more massive sources with stronger 
radiation fields. We modeled the source with the mea- 
sured fluxes in those bands, and also treating them as 
upper limits (since our model grid does not include PAH 
emission), and this difference did not greatly change the 
results. The source is not detected at J, so it must still be 
very embedded in its molecular cloud - indeed, the best- 
fitting models have Ay ~3G. Those models also have 
quite massive disks and envelopes, with accretion rates 
of a few times lO~^M0/yr. However, considering that 
this is an intermediate to massive star, that accretion 
rate does not necessarily imply extreme youth. Clearly, 
foUowup with spectroscopy and higher resolution pho- 
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tometry at longer wavelengths would be interesting to 
clarify the nature of this source, which may be respon- 
sible for some of the extra heating in that part of the 
molecular cloud seen in Figure [T] 




Fig. 20. — IRAC image of a massive young stellar object to the 
west of NGC 6611. Red, green, and blue are IRAC [3.6], [4.5], [8.0]. 



5. CONCLUSIONS 

We analyze Spitzer IRAC/GLIMPSE images of the 
M16 region. The diffuse emission in the long- wavelength 
IRAC bands (most likely PAH emission) traces the bub- 
ble blown by the young open cluster NGC 6611, and 
highlights the famous HST star-forming pillars. We use 
a large grid of model spectral energy distributions of 
young stellar objects to classify all the point sources in 
the area (<0.45° of Z=17.0, 5=0.8), and select candidate 
YSOs. The source classification determined using our 
model grid is comparable to classification based on IRAC 
colors (e.g Allen ct al. 2004; Mcgeath et al. 2004). Si- 
multaneous use of all photometric information between 
1 and 25 ^m provides additional insight into physical pa- 
rameters like envelope accretion rate and disk mass, and 
how tightly the data can constrain the circumstellar dust 
distribution. We particularly highlight sources that can 
only be fit well by models with massive envelopes, and 
are thus likely very young. The evolutionary classifica- 
tion possible using MIR data is not always the same as 
results based on NIR color excess alone. 

The distribution of young objects in the region is con- 
sistent with a picture of moderate distributed star for- 
mation, but no particularly vigorous activity triggered 
by the optical cluster NGC 6611. We confirm the youth- 
ful nature of that cluster, previously quantified as a 
large fraction of L- band (3.5 /im) excess sources (85% 
lOliveira et al."2005' ), and a large number of emission-line 
stars (52, ide Winter et al.iri997. ). Nearly all (>95%) of 
the sources in the NGC 6611 cluster that are detected in 
at least 4 of the 7 bands (J,H,K,[3.6],[4.5],[5.8],[8.0]) are 
YSO candidates. We identify a second cluster of YSOs 
to the northwest, apparently in a concentration of dense 
gas, and a cluster of YSOs at the tip of pillar 5 to the 
east. Examination of the two-point correlation function 
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Fig. 21. — Model fits to the western massive YSO. Plots and 
annotation as Fig |12l 

reveals that the YSO candidates are highly clustered. 
The nearest neighbor distribution suggests a character- 
istic size scale ~0.1pc, which corresponds to the Jeans 
mass in the clumps at the heads of the pillars. JVIid-IR 
imaging allows us to relate the YSO population to the 
physical conditions in their natal cloud and to paint a 
consistent picture of the star formation in JVI16. A newly 
discovered intermediate to massive YSO may be respon- 
sible for additional heating of the molecular cloud west 
of NGC 6611. 

Finally, we present a detailed examination of selected 
sources (for example in the HST pillars). Analysis of 
the 2-24 nui spectral index, bolometric luminosity, and 
all of the well-fitting model SEDs allows us to constrain 
the evolutionary states, stellar and disk masses, and ac- 
cretion rates, in some cases much more tightly than was 
possible from NIR data alone. Several sources in the 
trunks are very young/unevolved and it is possible that 
these sources were triggered by the compressive action of 
NGC 6611. The small ages of those sources could even 
be consistent with the very short destruc tion timescal e 
for the trunks by velocity gradients ( 1 0^ vr IPoundll 1 998h . 
On the other hand, some of the other sources in the 
trunks do not have massive envelopes; these may be more 
evolved, consistent with Class II or T-Tauri sources, or 
perhaps photoablation from NGC 6611 has accelerated 
the dispersal of their envelopes (both are located at pil- 
lar tips). If these sources are more than a few ~10^yr 
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old, it is difficult to understand how they could have 
been formed by triggering action of l-2Myr old stars in 
NGC 6611. We also do not see any particular pattern 
in evolutionary state, with younger sources farther from 
NGC6611, as might occur in a situation dominated by 
triggering (although we note that a pattern might be 
hidden b y the 3D str ucture of the pillars region as re- 
vealed bv lPoundlll998D . We favor the interpretation that 
YSOs formed spontaneously in dense parts of the cloud, 
and then these densest cores have persisted because they 
are more robust to photo-destruction than surrounding 
material; YSOs are found in exposed pillars, but were 
not necessarily formed as a result of that exposure. 
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TABLE 1 

Specific YSO candidates: observed parameters 



05 



name 


R.A. 


Dec. 


J 


H 




[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24] 


spectral index 


log{L.) 




(J2000) 


(J2000) 


[mjy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


2-24 ^tm 


log(L0) 


PI 


18hl6m00.1s 


-13d50mlls 


2.5±0.5 


6.6±1.0 


12. ±2.0 


26. ±4.0 


50. ±10. 


45. ±30. 


40. ±30. 


2000+1000 


1.0±0.2 


44±16 


P2 


18hl8m51.5s 


-13d49ml3s 


0.4±0.25 


1.0±0.3 


2.7±0.5 


4.3±0.5 


5.9±0.5 


10.6±2.0 


70±60 


<420=' 


0.3±0.2 


>5.3b 


Tl 


18hl5m58.5s 


-13d51m06s 


6.4±0.7 


12±1.2 


19±2.0 


28±8.3 


29±4.4 


43±8.6 


40±10 


560±400 


-0.3±0.1 


22±8 


P3 


18hl8m52.3s 


-13d49m38s 


1.0±0.7 


1.5±0.5 


2.2±0.5 


4.8±0.5 


17±2.0 


46±5.0 


70±40 


400±200 


2.3±0.2 


15±5 


T2 


18hl8m53.2s 


-13d49m57s 


2.8±0.5 


4.6±1.0 


5.9±1.0 


8.6±0.9 


12±1.2 


20±5.0 


40±30 


<300=' 


0.2±0.2 


6.0±1.5 


HH-N 


18hl8m58.9s 


-13d52m48s 


0.1±0.09 


0.3±0.2 


l.OitO.5 


3.0±0.5 


6.3±1.0 


6.0±4.0 


<4ia 


<310=' 


1.6±0.6 


>1.4b 


HH-S 


18hl8m59.4s 


-13d52m57s 


0.05±0.04 


0.2±0.1 


1.3±0.5 


1.7±0.5 


5.7±6.0 


6.0±2.0 


40±30 


<270=' 


1.1±0.5 


^3.0*^ 


P5 A 


18hl9m07.2s 


-13h45m24s 


0.6±0.3 


1.0±0.5 


4.0±2.0 


47±4.7 


170±44 


340±60 


440±80 


12000±10000 


2.0±0.1 


250±100 


P5 B 


18hl9m25.33s 


-13d45s35.5s 


0.5±0.3 


3.0±0.4 


24±2.4 


200±21 


320±150 


900±170 


1000±800 


2700±1000 


1.7±0.1 


200±50 


MYSO 


18hl8m08.6s 


-13d45m06s 


<0.2 


18±10 


140±100 


7000±500 


1200±800 


3700±1500 


5000±3000 


22000±10000 


1.0±0.3 


1000±300 



^ 3a" upper limit Without a measured flux density at 24 fim, YauF^ is certainly lower than the true bolometric luminosity. Since we do not include flux densities longward of 
MIPS-24 //m, all of the luminosities in this tabic may be low - this is reflected in the luminosities and masses of well-fitting models being somewhat (10-100%) higher. 



TABLE 2 

Specific YSO candidates: derived parameters 



name 


Av 


M* [Mq] 


log(M [Mo/yr]) 


log(Mdi,fc [M©]) 




min ave^ max 


estimate'' min ave max 


min ave max 


min ave max 



PI 

P2 
Tl 
P3 
T2 
P5A 
MYSO 



2.0=^ 
2.0 
2.0 
2.0 
2.0 
2.3 
27 



4.4 
14 
5.0 
5.9 
8.4 
6.8 
31 



7.0 

19 
5.7 
9.1 
11 
8.8 
33 



3 
1 

2.5 
1.5 
1.5 

6 



2.9 
0.4 
2.6 
0.9 
0.2 
0.5 
7.7 



4.6 
1.5 
3.2 
1.8 
3.2 
3.1 



6.5 
3.7 
5.3 
5.7 
5.5 
6.5 
9.8 



-4.8 



-5.3 



-4.9 



-4.5 
-5.2 
-6.1 
-4.3 
-6.3 
-4.5 
-6.2 



-4.2 
-3.5 
-4.8 
-3.2 
-3.7 
-4.3 
-3.8 



3.1 


-1.7 


-0.5 


4.0 


-2.5 


-1.3 


3.6 


-2.3 


-0.8 


4.0 


-2.9 


-1.2 


7.4 


-5.4 


-1.2 


2.6 


-1.6 


-0.7 


3.6 


-1.6 


-0.2 



p 
a- 

CD 

o 

% 



* For each parameter, the minimum and maximum value for any model which fits the data better 
than xl < 3.7 (probability of 0.001) is given, along with the average weighted by the probabihty 
exp(— x^/2).'' Mass estimated from an age consi stent with the 2-2 4 ^m spectral index, and the measured 
bolometric luminosity, using PMS tracks from [Sless et all II2000I) . This method is not expected to be 
more accurat e than a factor of 50% i n mass."^ Line-of-sight extinction was restricted to be greater than 
Av=2, since iHillenbrand et al.l l|1993f ) determined that the extinction to NGC6611 is Av'=3.1, and that 
cluster is in one of the less embedded parts of M16.'' Some disk-only models, with no accreting envelope, 
fit the data. 



